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High resolution N-body simulations have all but converged on a common empirical form for the 
shape of the density profiles of halos, but the full understanding of the underlying physics of halo 
formation has eluded them so far We investigate the formation and structure of dark matter halos 
using analytical and semi-analytical techniques. Our halos are formed via an extended secondary 
infall model (ESIM); they contain secondary perturbations and hence random tangential and ra- 
dial motions which affect the halo's evolution at it undergoes shell-crossing and virialization. 
Even though the density profiles of NFW and ESIM halos are different their phase-space density 
distributions are the same: p/a^ °^ r^", with a = 1.875 over ^ 3 decades in radius. We use two 
approaches to try to explain this "universal" slope: (1) The Jeans equation analysis yields many 
insights, however, does not answer why a = 1.875. (2) The secondary infall model of the 1960's 
and 1970's, augmented by "thermal motions" of particles does predict that halos should have 
a = 1.875. However, this relies on assumptions of spherical symmetry and slow accretion. While 
for ESIM halos these assumptions are justified, they most certainly break down for simulated 
halos which forms hierarchically. We speculate that our argument may apply to an "on-average" 
formation scenario of halos within merger-driven numerical simulations, and thereby explain why 
a = 1.875 for NFW halos. Thus, p/o^ °^ ^-1.875 ^ generic feature of violent relaxation. 
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1. Introduction 

There is a broad consensus that gravity driven evolution of the space distribution of collision- 
less dark matter, starting from some realistic matter power spectrum results in virialized objects 
whose spherically averaged density profile is well represented by the NFW prescription [^] or its 
close variants However, why the dark matter profiles have this shape is as yet to be determined. 
In an effort to shed light on the issue, [Q] investigated the phase-space structure of dark matter halos 
from N-body simulations. They found that the "poor man's" phase-space density, /(r) = p/a^*, is 
a power law over 3 decades in radius, /(r) oc r^", with exponent a = 1.875 = 15/8. 

We have developed an alternative method of generating dark matter halos, based on [^. It is an 
analytical scheme which treats collapse, shell-crossing and virialization of spherically symmetric 
halos. The halo contains secondary perturbations whose properties are calculated from the same 
power spectrum that gives rise to the main halo. These secondary perturbations induce random 
tangential and radial motions within the halo. The ensuing collapse can be likened to slow accretion 
of lumpy matter; there are no major mergers in our scheme; [^]. We call these halos ESIM, or 
Extended Secondary Infall Model halos. Even though ESIM and NFW halos are generated in very 
different ways, and have different density profiles, their phase-space density distribution, and the 
value of a is the same for both. Our goal is to understand why is a = 1.875 = 15/8 for both. 

2. The Jeans equation analysis of phase-space density distribution 

Equilibrium, non-rotating halos with isotropic velocity ellipsoids obey this Jeans equation: 

-r^ d{pd^)^ d 
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Following [Q] we use scaled variables, x = r/r^, and j = p/po- We assume power law behavior 
of /(r) with constant a for any given halo, but allow the density profile to have changing log-log 
slope, y oc jc^Pf'^). With these, eq. ^ can be rewritten in terms of x, y, exponents a and P, and a 
normalization constant. Unlike [0] we work with an equation obtained by differentiating eq. |2.1|: 



15p"-3p'(8a-5p-5) = (2a + p - 6)(2a- 5p)(2/3[a- p] + 1) (2.2) 
Here, p' and p" are derivatives of p with respect to — Injc. (Recall that p = —d\ny / d\nx.) 



Because eq. has many solutions (i.e. halo density profiles), it would help to make analytic 



inroads into the analysis of its solution space. We do just that: we derive a family of solutions: 

p' = p:„-2/3(p-po)^ p:„ = 3(2a + 5)2/200, po = (14a+ 15)/20. (2.3) 

All the halos that obey these three equations asymptote to very nearly constant P slopes at large and 
small radii, and have non-zero derivatives of P in-between. We will refer to this as the main family 
of solutions. Each member is completely defined by specifying a. As far as we can tell, this is the 
only family that has a closed form analytical description. Note that the 'critical solution' identified 
by [0] is a member of this family; it is obtained by setting a = 1.875. 

Figure |l|(a) is a schematic representation of all solutions of eq. 2.2. A few main family so- 



lutions are shown as solid vertical red line segments; the inner density cusp for each is given by 
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Figure 1: (a) Pictorial depiction of solutions of eq. 2.2. a — 1.875 is the slope of the phase-space density 
profiles of NFW and ESIM halos; a = 1.9444 results in periodic, zero-damping solutions of eq. 2.2. Six 
solid vertical line segments are examples of the main family of solutions. Labeled lines: factors in eq. 2.2; 
Po is given in eq. 2.3. (b) A set of solutions of eq. 2.2, with constant a = 1 .875 but varying initial values of 
P = 1.25, 1.75, ...3.25; initial p' = for all. The thick line belongs to the main family of solutions; note that 
NFW empirical fit follows this solution closely for intermediate values of p. Dotted and dashed lines are two 
ESIM halos; short vertical line segments indicate their virial radii. Normalization of radius x is arbitrary. 



(2a — 5P,„) = 0, the outer density profile slope is given by (2/3 [a — ^out] + 1) =0. Fig. |T](a) illus- 
trates the role these factors on the RHS of eq. ^!2|play in defining main family solutions. 

One of these six vertical line segments is for a = 1.875; the corresponding P(x) is shown as 
thick solid curve in fig. |l](b). Other curves are what we call a-family solutions, obtained by keeping 
a constant, but varying initial conditions: P and P'. Note that one of the members of the a-family 
is a power law density profile, P=const. The P value for the power law solution of an arbitrary 
a-family is obtained by solving (2a + P — 6) = 0; see eq. 'LI and long-dash line in fig. |l|(a). 

Does eq. 2^ allow any special a values? The equation is a non-linear damped oscillator. For 
some a the x-averaged value of the 'dissipation' term, — 3P'(8a — 5P — 5) is zero. This value is 
a = 35/18 Rr: 1.9444. When a = 35/18 any initial condition for P and P' results in a periodic, 
constant amplitude function of p vs. x. Value a = 35/15 is only ~ 4% different from a = 1.875. 
Aside from this coincidence (or not?), and many other interesting insights provided by the analysis 
of eq. [2.2[ , we were not able to identify why a = 1.875 is special. Next, we take a different, and 
more promising approach to uncover why a = 1.875 for ESIM and NFW halos alike. 



3. Secondary infall analysis of phase-space density distribution 

Consider the initial stages of halo collapse in the context of secondary infall models. A small 
constant central mass excess, 6Mo is surrounded by material of average density. The dynamics of 
the pre-tum-around period is described by the parametric equations of The turn-around radius 
for a shell of initial comoving r,, is r„, = r,5^\ and 5 is the initial fractional overdensity inside r,. 
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Then, M(< r,) oc pr?, so r,- oc M^/^ . Also, 6 = 5Mo/M oc M^^. Combining these scalings we get 
r,„ oc M"*/^. After reaching turn-around each shell collapses back a little; one typically assumes a 
constant collapse factor. Assume that each shell spends most of the time at its apocenter, and so 
most of its mass is located at that radius, which is oc r,„. The resulting density distribution in the 
proto-halo is, p(r) oc r^^dM/dr, or, p oc M^^ oc rffi . This is a well known result. 

In the real Universe the collapse of material will not be purely radial; there will be some 
random motion of particles, and associated kinetic energy. We speculate that during the early stages 
of collapse the kinetic energy will be derived from the gravitational potential energy, and therefore 
will be proportional to the potential energy: ^v^ oc GM/r,„. Kinetic energy of random motion gives 
an estimate of the velocity dispersion: a oc [M/r^]''^ ' . So, . Combining p(r) 

from the previous paragraph with o{r), we get p/o^ oc (r^^/4)/(r^3/^) oc r^^^/^^ i_e. a = 1.875! 

The result obtained above says that the phase-space density is a function of r. Because 
E oc GM/r, £" is a monotonic function of r, at least for systems that are, on average, spherically 
symmetric. In equilibrium, the total energy £ of a particle is its integral of motion. If the halo 
collapse proceeds slowly then the halo passes through a series of quasi-equilibrium stages. Maybe 
we can assume that energy is an 'approximate' integral of motion in a slowly collapsing galaxy. 
In that case, our relation p/a^ oc ^- 1-875 interpreted as saying that the phase-space density 

(p/a^) is a function of E only, in compliance with Jeans Theorem. 

If dark matter is collisionless, then the collapse will preserve the phase-space density as calcu- 
lated above. So the final virialized dark matter halos, whatever their density and velocity dispersion 
profiles, will have the same phase-space density distribution that was characteristic of the early 
stages of collapse. Therefore, virialized halos are expected to have p/a^ oc 875 jj^^ argument 
relies on many approximations, most of which can be justified for spherically symmetric, smoothly 
accreting ESIM halos. However, one is hard pressed to see why these approximations would hold 
in a hierarchical formation model. We speculate that the argument could apply to an 'average' 
situation; after all, NFW is an average shape of numerically generated halos. 

Finally, we note the connection between the above argument and the a-family solutions of § 0. 
The various evolutionary stages of halos have different density profiles. We claim that all these are 
represented by members of the same a = 1.875-family. The earliest epoch halos have power law 
density profiles with P = 9/4, value derived at the top of this §; the same value is derived using 
Jeans equation analysis — the horizontal line in fig. [I](b). Later stages of halo evolution develop 
changing slope density profiles, represented by more and more curved lines in fig. [I](b). The final 
stage, the thick solid line is a good approximation to the NFW profile obtained from simulations. 
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